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Twist-3 distribution amplitudes of the pion and kaon are studied in this paper. We 
calculate the fist several moments for the twist-3 distribution amplitudes (</>p jCT and 
4>p,a) °f the pion and kaon by applying the QCD sum rules. Our results show that, (i) 
the first three moments of 4>p and the first two moments of and 4>Z' K of the pion 
and kaon can be obtained with 30% uncertainty; (ii) the fourth moment of the <j>p 
and the second moment of the can be obtained when the uncertainty are relaxed 
to 35%; (iii) the fourth moment of the can be obtained only when the uncertainty 
are relaxed to 40%; (iv) we have = 1.10 ±0.08 GeV and m\ K = 1.25 ±0.15 GeV 
after including the a; s -corrections to the perturbative part. These moments will be 
helpful for constructing the twist-3 wave functions of the pion and kaon. 
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I. INTRODUCTION 



Hadronic distribution amplitudes, which involve the non-perturbative information, are 
the important ingredients when applying the QCD to hard exclusive processes via the fac- 
torization theorem. These distribution amplitudes are process-independent and should be 
determined by the hadronic dynamics. They satisfy the renormalization group equation and 
have the asymptotic solutions as Q 2 — > oo. 

From the counting rule, the twist-2 distribution amplitude makes the leading contribution 
and the contribution from the higher-twist distribution amplitude is suppressed by a factor 
1/Q 2 in the large momentum transfer regions. However as one wants to explain the present 
experimental data, the non- leading contributions should be taken into account. The non- 
leading contributions include higher-order corrections, higher-twist and higher Fork state 
contributions et.al. Therefore one has to study the twist-2 and higher-twist distribution 
amplitudes as the universal nonperturbative inputs for the exclusive processes. 

Distribution amplitudes can be obtained from the hadronic wave functions by integrating 
the transverse momenta of the quarks in the hadrons. For example, the pionic distribution 
amplitudes of the lowest Fork state are defined as : 

(0 \d a (z) [z, -z\ up(-z)\ ir(q)) 

= ~ g/ w J i # | foMZ) + <r750£(O + jjm^ 7 5g^€(0}^ + ■■■(!) 

where a^ u = 1(7^7*, — 7^7^), fw is the pion decay constant and 

[z, —z] = exp sig dx^A^ 



is the Wilson line inserted to preserve gauge invariance of the distribution amplitudes. The 
071- (0 , 4>p(0 an d 4>a(0 m Eq.((TJ) are the twist-2 and two twist-3 (non- leading) distribution 
amplitudes respectively. For the K meson, the definition is similar except for the d quark 
replaced by the s quark and replaced by . 

To isolate the light-cone twist-3 distribution amplitudes 0£ and 0^ of the pion, one can 
contract the equation (JTJ) with the gamma matrices 75 and 07^75 respectively, 

<0 \d{z)i lb [z, - Z ] u(-z) I n+(q)) = m p J« \ J d£ + '~ ( 2 ) 



and 

<0 |d»<V75 [z, -z] u{-z)\n + {q)) = ~ l m ^ h {q,z u - q v z») \ £ dt, <%{t)e«M 

+ ■■■■ 1 (3) 

In a similar way, we can define two twist-3 distribution amplitudes (ftp and ft of the kaon 
in the following, 

<0 \s(z)i l5 [z, -z] u{-z)\K + {q)) = m p QK f K I J d( <%(Qe*™ + ■■■ (4) 



and 



{Q\s{z)v^[z-z}u{-z)\K + {q)) = f^{q,z v - q„z,) -J d( ft (Qe^ 

+ ■■■■ ^ (5) 

The dots in the above definitions are refer to those higher-twist distribution amplitudes. We 
do not consider their influences in the following calculation. 

As pioneered in Ref. , the authors pointed out that the first several moments of the 
distribution amplitudes could be calculated in the QCD sum rules Those moments are 
helpful to construct a model for the hadronic distribution amplitudes. 

The parameters and vnP^ introduced in the definition are used to normalize the 
zeroth moments of their corresponding distribution amplitudes. It is shown in this paper 
that these parameters determined by the QCD sum rules are smaller than those required by 
the equations of motion (e.g., see Ref. 0,0])- 

In this paper, we calculate the first three moments of twist-3 distribution amplitudes of 
the Ti and K , defined in Eqs.©-©, in the QCD sum rules. For the pion case, we had 
calculated the moments of distribution amplitude 0£ in the previous paper However 
there were some mistakes in estimating the contribution from the continuous spectrum 
and the Borel windows which would severely influence the values of the moments. Now 
we present the correct expressions for the moments of 0£ and re-analyse their numerical 
results in this paper. Furthermore, it is well known that axial currents in a correlator would 
couple to instantons(see, for example, Ref.|7|). And this may cause some complication in 
the calculation and make the results unreliable. We will not explore their influences in this 
paper. 



This paper is organized as follows. In Sec. II, we give the sum rules of the moments of 
4>p and <f)^ for the 7r meson. The sum rules for the moments of (f)^ and 4>^ of K meson 
are given in Sec. III. The SU(3) symmetry violation have be taken into account. In Sec. IV, 
numerical analysis of various moments is presented. The information about the 3-particle 
twist-3 distribution amplitude obtained from the 2-particle distribution amplitudes are also 
discussed. The last section is reserved for summary and discussion. 



II. QCD SUM RULES FOR THE MOMENTS OF 0£ AND #J OF THE PION 



In this section we apply the background field method in QCD to calculate the moments 
0. floL llll ll^ . Expanding equations (j2J) and (jHJ) around z 2 = 0, we have 



<J(0)7b(« • D )»u(0) n + (q)\ = -if^ <^> (z ■ q) n 



(6) 



and 







d(Q)a^ 5 (iz- D ) n+1 u(0) n + (q) 



n + 1 



Um^(C)(q^-q^)(z-qr (7) 



respectively. The moments in the equations (jBJ) and (JJJ) are defined by the following expres- 
sions, 



(8) 



As usual, the SU(2) isospin symmetry can be taken as (nearly) exact. It means that the 
distribution of longitudinal momentum carried by the quarks (in the light cone framework) 
should be symmetric between u and d, i.e., odd moments of the distribution amplitudes 
<pp, (ftp should be zero. So we consider only the even moments for the pion case in the 
following. 

In order to obtain the sum rules of the moments, we introduce two corresponding corre- 
lation functions, 



(,.^/M ((?2) 



d A xe iqx ( 



T\d{x) l5 {iz- Dy n u(x),u(0)-f 5 d(0) 







(9) 



and 



-i(q,z u - q u z,)(z ■ q) 2n I^°\q 2 ) 
-/ / r/'./v'"-' /() 



T [d(x)a^j B (iz ■ D Y n+1 u(x), w(0) 75 d(0) 



(10) 



In the deep Euclidean region( — q 2 3> ), one can calculate the Wilson coefficients in the 
operator product expansion (OPE) for Eqs.(jn} and (fTUj) perturbatively. The results with 
power correction to dimension six and the a s corrections to lowest order are written as 



i {2nfi \ q 2 



1 



)QCD : 
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3 2l -q 2 
q in 
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—G 2 



2n - 1 (m u + m d ) (qq) 



and 



Il 2n '°\q 2 



)QCD : 
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3 2] ~q 2 



8 q 2 

® s qq) 2 

q' 



(11) 



a 



S G 2 



7T 



2n + 1 (m u + m d ) (gg) 



2n + 3 8tt 2 



167T , 9 

+ — (16n 2 + Yin - 7) 



24 q 2 

2 



__. . -. (12) 

81 v g 4 

On the other hand, in the physical region, correlation functions (JOJ) and (jlOj) can be written 
in terms of their hadronic spectrum representation (according to Eqs.fJOJ) and (fTT|). Q and 
(|T2*|) respectively), 



Im I^°\q\ ad = vr5(g 2 - m 2 )/ 2 (m^) 2 <£f ) + vr- 



\ix 2 2n + 1 



g 2 % 



2 



(13) 



and 



Im I^°\q\, d = n5(q 2 - m 2 f-^f 2 m^ + vr- 3 



-q 2 9(q 2 - O . (14) 



The correlation function in these two regions can be related by the dispersion relation, 

1 



7T 



qcd 



In order to improve its convergence, we apply the Borel transformation, 

\w I ds e ~ S/M2 Im /(s)had = Lm /( ~ g2) Q CD ' (15) 

where M is Borel parameter. Substituting (|TT|) and (fTB^) into (fT3j) gives the sum rules for 



the moments of 0^ : 



P \2. 
On) " 
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(2n + 1) 8tt 2 
2n - 1 (ra„ + m d ) {qq) 



8 M 4 
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(16) 



Similarly, substituting Eqs.flHJ) and (fT3j) into Eq.(fT5)l gives the sum rules for the moments 
of 61 : 
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(17) 



where s£ and are the threshold values to be chosen properly, and the zeroth moment 
has been normalized to unit , i.e., = = 1 . 



III. QCD SUM RULES FOR THE MOMENTS OF 



AND 0f OF THE KAON 



For the kaon, we should consider the difference between s quark and u quark (i.e., the 
violation of the SU(3) flavor symmetry). There is an asymmetry of the distribution of the 
longitudinal momentum carried by s quark and u quark in the light cone framework. So the 
odd moments of distribution amplitudes for K meson do not vanish. The violation effects of 
the SU(3) flavor symmetry for leading-twist distribution amplitudes of K and/or K* meson 
were considered in Ref. [14] . So in calculating the odd moments, we retain all the corrections 
to order m 2 . 



Expanding equations (JH) and (0) around z 2 = 0, one obtains 



5(0) 75 (w " D ) n u(0) K + (q)) = -if K m p QK (Q) (z ■ q) n 



(18) 



and 



s(0)<V75(is- D) n+1 u(0) K + (q) 



n + 1 



fKm° K (C)(q,z u -q u z,)(z-qr (19) 



respectively, and the moments are defined by 



(20) 



2 y_i Tp w ' ' ^ a< 2 j _ j 

Similar to the pion case, the correlation functions for calculating the moments of the kaon 
are defined as, 



[z-q) n I { ^\q 2 ) = -i I d 4 xe^ x (0 



T s{x) l5 {iz ■ D ) n u(x), u(0) 75 s(0) 







(21) 



and 



-i{q^z v - q v z^){z ■ q) n p£ g °\q 2 ) 



d 4 xe iq x ( 



T {s(aO<v7 B (is ■ *D ) n+1 «(x), ^(0) 75 s(0)} 







(22) 



As discussed in the previous section, the correlation functions can be calculated perturba- 
tively in deep Euclidean region, i.e., Q 2 = —q 2 ^> 0. Combined with Eqs. (|18j) and ()19jl . we 
assume the hadronic spectrum representations of the above correlations as follows, 



1 



and 



;<? 2 % 2 



(23) 



(24) 



Employing the dispersion relation and Borel transformation as done in previous section, the 
sum rules for the moments of can be expressed in the following, 
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(25) 



and 



2 /M 2 M 4 



/ 2 



A" 



3 m 2 
"8tt 2 M 2 
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(26) 



For the twist-3 amplitude 



we make a similar calculation according to the above 



procedure and the sum rules for the moments of <ft^ become 
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(27) 



and 



r i\ m a _3 e^/ M2 M 4 f 1 m 2 . / 4/ ^\ , m« (ss) - m u (uu) 



M 4 



6M^ (ln ^ + 1 ~ 7£)_ 3 
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27 M 6 



(28) 



where 7^ = 0.577 216 • • • is the Euler constant, s p K and in the above equations are 
the threshold values to be chosen properly, and the zeroth moments has been normalized, 

<C P °> = (C) = 1 • 



IV. NUMERICAL ANALYSIS 

To analyse the sum rules |T6|l . |T7j) . (j25|) - (j28|) numerically, we take the input parame- 
ters as usual : f K = 0.160 GeV, f n = 0.133 GeV, m s = 0.156 GeV, m u = 0.005 GeV, 
m d = 0.008 GeV, (uu) = (dd) = -(0.24 GeV) 3 , (ss) = 0.8 (uu) , g s (saGs) = 
-0.00885 GeV 5 , (—GG^ = 0.012 GeV 4 , a s (lGeV) = 0.5 . The renormalization scale 
/i = M is assumed in the following analysis. 

As to the threshold values s v ^ a K in the sum rules, they can be taken to the mass square of 
the first exciting states in the corresponding channel. Although the windows become broader 
when the s p ^ u K are larger, the threshold values can not exceed the first exciting states. In 
order to get maximum stability of the sum rules, they are taken as the mass square of the 
first exciting states, i.e., 

= (1.3 GeV) 2 , s p / = (1.46 GeV) 2 

where the first exciting state is 7r'(1300) for the pion case, and that for kaon case is K(1460) 



A. Determination of the normalization constants 



For each distribution amplitude, we introduce a corresponding parameter (e.g., for 
4>p). These parameters are normalization constants which normalize the zeroth moments to 



one. Their values can be determined from the sum rules (|16j) . (|17jh (|25|) and (|27|1 with n = 0. 
Take m^. as an example. To identify a Borel window (M 2 ) for the sum rule of m^, one 
requires that the continuum contribution is less than 30% and the dimension-six condensate 
contribution is less than 10%. This requirement leads to a window M 2 E (0.64 , 0.75) GeV 2 
and one can find = 0.96 ± 0.03 GeV within this window. And the results are plotted in 
the Fig da) and Fig Kb) for 
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FIG. 1: (a) The window for the normalization constant m^L without a s correction in the perturba- 
tive part in the sum rule. The dashed line is the ratio of the dimension-six condensate contribution 
to the total sum rule(n = 0) and the solid line is the ratio of the continuum contribution to the 
total sum rule(n = 0) . (b) The corresponding values of rriQ n within the window. 



The same procedure can be applied to get other parameters m^, mj^, the numerical 
results are listed in Table |l] . The continuum contribution to the sum rules are required 
to be less than 30% and the dimension-six contribution is required to be less than 16% for 
m^, 10% for m p 0K and TOq K . It should be pointed out that when the «, correction to the 
perturbative part of the sum rule for m p QlTl fn p QK are taken into account |9J, their values will 
be increased by 15-20% . For example, m p 0n = 1.10 ± 0.08 GeV and m p 0K = 1.25 ± 0.15 GeV. 

One can see from above that is smaller than about 30%. The main reason is 
that the opposite sign of the dimension-six condensate terms in Eqs. (II 6|) and (|T7j) . For 
the kaon case, the approximate 30% difference between vn^ K and vnP QK is due to the same 
reason(see Eqs. (j23j) and (jUj)). 

It was shown that the normalization constants for the twist-3 distribution amplitudes 
can be obtained from equations of motion 4]. So at this point, we would like to compare 
our results with those obtained by equations of motion and judge upon the accuracy of 
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m (GeV) 


0.96 ± 0.03 


0.67 ±0.06 


1.06 ± 0.09 


0.71 ±0.09 


M 2 (GeV 2 ) 


0.64-0.75 


0.60-0.68 


0.58-0.93 


0.66-0.83 



TABLE I: Normalization constants itiqS and the corresponding Borel windows for the distribution 
amplitudes <fjl a and without a s correction in the perturbative parts in the sum rules. 

the sum rules presented above. From the equations of motion follows the normalization 
constant — > fi n = m 2 / fn for (p^ and — > Jl % = [i n — m for 0£ , where m = m u + 
m,i. Furthermore, the Q 2 dependence of the quark mass fh can be written as: rh(Q 2 ) = 
[ln( / u 2 /A 2 )/ln(Q 2 /A 2 )] 4//9 m(/i 2 ) which can be obtained from the anomalous dimension of 
quark mass. Thus we find fj, n (l GeV 2 ) ~ 1.48 GeV and GeV 2 ) ~ 1.47 GeV as we take 
m(4GeV 2 ) = HMeV. For the kaon case, we take (m u + m s )(4GeV 2 ) = 140 MeV. From 
the equation of motion we have /ix(lGeV 2 ) = m 2 K /(m u + m s )(lGeV 2 ) ~ 1.45 GeV and 
(/i^-(m s ±m u ))(lGeV 2 ) w 1.28 GeV. In the above statement, QCD scale A = 250 MeV is 
assumed and n/ = 3 flavors are taken into account. One can see that the deviation of 
from jjtTx is aboat 26% and the deviation of m^ K from hk is less than 15% (a s corrections to 
the perturbative parts are included). 

If the a s correction to the perturbative part is 15-20% and these corrections make the 
normalization constants increasing, one expect that the deviation of from Jl n is about 
45% and the deviation of mjj^ from \ik — {ra u + m s ) is about 33% . 

B. Determination of the second moment of <$^ a and the odd moment of <pp a 

Let us consider the second moments of cf) 1 ^ and <P1 for the pion. Just as the determi- 
nation of normalization constants in the above paragraphs, one should find a window for 
each moment in the corresponding sum rule. The Borel windows in Table |H] are obtained 
under requirement that both the contributions from continuous states and the dimension-six 
condensate are less than 30% . As an example, we plot the results for the moment (£ 2 ) in 
the Fig 12(a) and the Fig 12(b) and the numerical results are listed in Table |H] 

Now we turn to the determination of the first moments of 4>p a - The contribution from 
the dimension-6 condensate and the continuous states of (Q) and are plotted in FigH3 
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(a) (b) 
FIG. 2: (a) the window for the moment the dashed and the solid line indicate the ratio of the 
contribution of dimension-6 condensates and continuous states in the total sum rule respectively; 
(b) the moment within the Borel window. 

For , the dimension-six contribution is less than 1% and the continuum contribution 
is less than 10% . For (Q) , the contribution of dimension-six condensate and continuous 
states are less than 10% . 

With these windows we can get the values of the corresponding moments. These results 
are listed in Table HH 
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FIG. 3: The windows for and (Cs-)- The dashed and the solid lines indicate the ratios of the 
contribution of dimension-6 condensates and the continuous states in the corresponding total sum 
rule respectively. 
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0.52 ±0.03 


0.34 ± 0.03 


-0.10 ±0.03 


-0.13 ±0.04 


M 2 (GeV 2 ) 


0.72-0.88 


0.71-0.84 


0.80-1.85 


0.77-1.53 



TABLE II: Second moments of odd moments of (pp a and their corresponding Borel windows. 
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0.44 ± 0.01 


0.20 ±0.01 


0.43 ± 0.04 


0.173 ±0.002 


M 2 (GeV 2 ) 


1.06-1.14 


1.08-1.22 
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TABLE III: Fourth moments of (j)l (J and second moments of 4>p a - But note that the values of 
(£p)> (Cct) given in this table is under the requirement of 35% uncertainty and is under the 
requirement of 40% uncertainty. 

C. Determination of the fourth moment of 0£ and the second moment of <pp a 

Now we consider the second moment (£ 2 ) of 4>p for the K meson. The Borel window for 
(Cp) is shown in the Fig0Ja) when the contribution of continuous states and the dimension- 
six condensate are less than 30% . The numerical results are listed in Table IIHI 

However, for the fourth (n = 2) moments of <j)^ of the n meson, we can not find the 
Borel windows when the contribution of continuous states and the dimension-six condensate 
are required to be less than 30% . For the second (n = 1) moment of <fi^ of the K meson, we 
find that the Borel window is very narrow when the contribution of continuous states and 
the dimension-six condensate are required to be less than 30% . As we relax the requirement 
that the contribution of continuous states and the dimension-six condensate are less than 
35%, the Borel windows for (£ 2 ) and can be found. For , one can find the Borel 
window only when the contribution of continuous states and the dimension-six condensate 
are less than 40% . The above windows are shown in Fig. E^b)-0Jd) . 

The values of these moments within their corresponding windows are listed in Table IIHI 
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dashed and the solid lines indicate the ratios of the contribution of dimension-6 condensates and 
the continuous states in the corresponding total sum rule respectively. 

D. From two-particle distribution amplitudes to three-particle distribution 

amplitude 

There are three twist-3 distribution amplitudes 0£ , <ft^ and (f)^ for tc meson. As shown 
in Ref.JJ] , they are not independent. By employing equations of motion in QCD, one 
can obtain some relations between them. For the pion, the relations between the twist-3 
distribution amplitudes of the lowest fork state and the 3-particle one are given in Ref-j^] . 
They obtained two distribution amplitudes of the lowest fork state <j)p a from the 3-particle 
distribution amplitude 4>3 7T which was given by a direct calculation in the QCD sum rule 
method . Contrarily, we use the relations from equations of motion to see what we can 
say about the 3-particle distribution amplitudes with the above results of the distribution 
amplitudes of the lowest fork state as input. The results can also be compared with those 
obtained by QCD sum rule directly[3] . The cross checks in these calculations are helpful 



to judge upon the accuracy of the sum rules. 

First, let's discuss the pion case. The three-particle distribution amplitude of the ti meson 
can be defined as, 

(0\d(x)a^ u 'j 5 gG a f 3 (-vx)u(-x)\7c + (q)) 

= ifz* [fo Ww* - gj^) - (a <- (3)] J $> ai e^-^+^^i^) (29) 

where 3>ai = daida2da^8{a\ + a 2 + «3 — 1) • There are the system of recurrence relations 
for the moments and : 
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n -\- o n + o j_i 

-6^ -/ d«v<((ar 2 -ai + va 3 ) n ~ 1 )) (31) 

where (((«2 — c*i +^«3) n )) = / 03 7r ( Q; i) {oi2 — Oii + va 3 ) n defines the moments of 3-particle 
distribution amplitude. Instead of taking R p = R a = R as in Ref. j^J , we introduce them 
seperately, 

-"-P — — p~~ ~7~ ancl -n-cj — — ~~F~ 
m 0n Jir m 0ir Jtt 

Now, we take second moments into account. The above relation can be reduced to 

<a = ~<a+f^-~iu«3>> 

and 

which gives, from Table H] and Table ITT1 

((a 3 )) = (0.13 , 0.27) and f 3n = (0.0049 , 0.0067)GeV 2 . (32) 
At this point, we compare the moment ((a^)) and the with those calculated directly by 



the sum rule method in Ref. 



lfj : ((a 3 )) = (0.06 , 0.22) , / 37r « 0.0035 GeV 2 . One can see 



that the results from the two approaches are compatible to the order of magnitude. 



From the analysis in previous section, we have shown that the fourth moments and 
can not be obtained in a reliable way, so we do not use them to give the other moments, 
i.e., {{af}) and ((aia^)) , etc.. 

Now we turn to the K meson case. Similar to the pionic case, one can define a three- 
particle distribution amplitude 

(0\s(x)a fM/ 'y 5 gG a p(-vx)u(-x)\K + (q)) 
= ifsK[q«(q?Sv{} ~ qJtf) ~ (a <- P)\ [ ®o« e^-^+^Wte) ■ (33) 



Following Ref.Q, a similar recurrence relation can be obtained. As the first and second 
equations, 



moment of are taken into account, the recurrence relation can be truncated to three 



(O = (Q (34) 

(£) = | f (O - l 5 K((">)U (35) 

«J> = (36) 

where i^, = fzx / {fK^xix) > an d the primes on the i?s and the subscript K indicate that 
the quantities are related to K meson. From Table HI we have R' a /B! p w 1.06/0.71 , so the 
equation (|H4jl is a direct constraint of the two fist moments. Our calculation (see Table HTjl 
shows that the left hand side of Eq. (|H4|) is about —0.13 and the right hand side is about 
—0.11. It can be seen that this equation is approximately fufilled. Solving the last two 
equations and (jHHj) , we can obtain f 3K and ((0:3)) k'- 

f 3K = (0.0071 , 0.0105) GeV 2 , ({a 3 )) K = (5.01 , 5.37) (37) 

To determine more moments of the three-particle distribution amplitude, we have to include 
higher moments of the two-particle distribution amplitudes. However, one can not guarantee 
the convergence of the operator expansion for bigger n . 

V. SUMMARY AND DISCUSSION 

In this paper we calculate the first three moments of the twist-3 distribution amplitudes 



and (fit; a by using the QCD sum rules. It has been shown that the first three moments 



of (ftp and the first two moments of cft 1 ^ and <ft^ ,K of the pion and kaon can be obtained with 
30% uncertainty. The fourth moments (£p jCr ) of (ftp <T and the second moment of (ft^ 
can be obtained under 35%-40% uncertainty. When the a s corrections (we take them from 
Ref. ^|) to the perturbative part of m p 0K , are included, we find that the values of 
and TTT-q^- £1X6 increased (and the corresponding Borel windows becomes a little narrower) 
to m>Q K = 1.25 ± 0.15 GeV and mf^ = 1.10 ± 0.08 GeV. It may be expected that the a s 
corrections to the perturbative parts in the sum rules for m% K and will be about 15-20%. 

As to the normalization constants m p Q '" and vnF^ , our calculated results show that they 
are smaller than the values which are given by the equations of motion and at the same 
time, the calculated m^. K are smaller than the corresponding K . These deviation can 
be traced to the non-perturbative condensate effects (see the sum rules (1TT)|) . (fTTj) . (|23j) and 
(|2Tj) for the normalization constants), in particular, the dimension-six condensate terms in 
opposite sign lead to about 30% difference between these normalization constants. On the 
other hand, from the sum rules, one can see that the contribution from the continuous 
state grow too fast which prevent us from taking M 2 to be larger values(larger M 2 will 
lead to bigger values of the normalization constants), and then the window for the sigma- 
sum-rules(mQ 7r K ) are much narrower than the non-sigma sum rules {m p Q7t K ). So we think 
the smaller values of K may be related to our approximation in the hadronic spectrum 
representation. 

Furthermore, we calculate the moments of the quark-antiquark-gluon distribution ampli- 
tude from the numerical results on the distribution amplitudes of the lowest Fork state b 
applying the exact equations of motion and compare our results with those from Ref. 
The comparison shows that they are compatible with each other to the order of magnitude. 
It is helpful to improve the accuracy of the QCD sum rules approach for getting more precise 
information on the twist-3 distribution amplitudes. 

These moments can provide several constraints upon the twist-3 distribution amplitudes. 
These constraints will be helpful for building the model of the distribution amplitude. For 
example, Ref suggests a model for the twist-3 wave function of the pion based on the 
QCD sum rule calculation to get a more realistic contribution to the pion form factor. Here 
we discuss the distribution amplitude (ft^ of the kaon(since the first three moments can 
be obtained reliablely). As usual, we expand the distribution amplitudes in Gegenbauer's 
polynomials and use the moments to determine their first few coefficients in a truncated 



; e by 



form : <j>p(C) = En=o^ 2 (O a n- From the three moments of (see Table ITT1 and Table 
we have the twist-3 distribution amplitude approximately, 

0f(C) = 1 - 0.30 C X 1/2 (C) + 0.73 C 2 1/2 (C) (38) 
which is asymmetric since the SU(3)f symmetry broken are taken into account. 

Acknowledgments 

We would like to thank A. Khodjamirian for his helpful communication. This work was 
supported in part by the National Science Foundation of China, 10275070. 



[1] V. L. Chernyak and A. R. Zhitnitsky, Nucl. Phys. B201, 492 (1982). 
[2] V. L. Chernyak and A. R. Zhitnitsky, Phys. Rep. 112, 173 (1984). 

[3] M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, Nucl. Phys. B147, 385 (1979); M. A. 

Shifman, A. I. Vainshtein and V. I. Zakharov, Nucl. Phys. B147, 448 (1979). 
[4] V. M. Braun and I. B. Filyanov, Z. Phys. C 48,239 (1990) 239. 
[5] P. Ball, V. M. Braun, Y. Koike, K. Tanaka, Nucl. Phys. B529,323 (1998). 
[6] T. Huang, X. H. Wu, M. Z. Zhou, Phys. Rev. D 70, 014013 (2004). 
[7] L. Baulieu, J. Ellis, M.K. Gaillard and W.J. Zakrzewski, Phys. Lett. 77B, 290 (1978); 

B.V. Geshkenbein and B.L. Ioffe, Nucl. Phys. B166, 340 (1980); 

V.A. Novikov, M.A. Shifman, A.I. Vainshtein and V.I. Zakharov, Nucl. Phys. B174, 378 
(1980); 

M.S. Dubovikov and A.V. Smilga, Nucl. Phys. B185, 109 (1981); 
B.L. Ioffe and A.V. Smilga, Phys. Lett. 114B, 353 (1982); 
V.A. Nesterenko and A.V. Radyushkin, Phys. Lett.ll5B, 410 (1982); 
E.V. Shuryak, Nucl. Phys. B203, 116 (1982). 
[8] V. A. Novikov, M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, Fortschr. Phys. 32, 585 
(1984). 

[9] L. J. Reinders, H. R. Rubinstein and S. Yazaki, Phys. Rep. 127, 1 (1985). 
[10] X. D. Xiang, X. N. Wang, T. Huang, Commun. in Theor. Phys. 6, 117 (1986). 



[11] T. Huang, X. N. Wang and X. D. Xiang, Phys. Rev. D 35, 1013 (1987). 
[12] V. Elias, T. G. Steele and M. D. Scadron, Phys. Rev. D 38, 1584 (1988). 
[13] T. Huang and Z. Huang, Phys. Rev. D 39, 1213 (1989). 

[14] A. Khodjamirian, Th. Mannel and M. Melcher, Phys. Rev. D 70, 094002 (2004), 
hep-ph/0407226| ; 

V.M. Braun and A. Lenz, Phys. Rev. D 70, 074020 (2004), |hep-ph/0407282] . 
[15] Particle Data Group, Phys. Rev. D 66, 585 (2002). 

[16] A.R. Zhitnitsky, I.R. Zhitnitsky, V.L. Chernyak: Yad. Fiz. 41, 445 (1985). 
[17] T. Huang, X. G. Wu, Phys. Rev. D 70, 093013 (2004). 
[18] P. Ball, J. High Energy Phys. 01, 010 (1999). 



